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Recent findings of superconductors that simultaneously exhibit multiple spontaneously broken symmetries, 
such as ferromagnetic order or lack of an inversion center and even combinations of such broken symmetries, 
have led to much theoretical and experimental research. We consider quantum transport in a junction consisting 
of a ferromagnetic metal and a non-unitary ferromagnetic superconductor. It is shown that the conductance 
spectra provides detailed information about the superconducting gaps, and is thus helpful in determining the 
pairing symmetry of the Cooper pairs in ferromagnetic superconductor. 
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I. INTRODUCTION 



Spontaneous symmetry breaking in condensed matter sys- 
tems ranks among the most profound emergent phenomena in 
many-body physics. Multiple spontaneously broken symme- 
tries are not only of interest in terms of studying properties 
of specific condensed matter systems, but also due to the fact 
that it may provide clues for what could be expected in other 
systems in vastly different areas of physics. A first attempt 1 at 
describing the physics of magnetic spin-singlet superconduc- 
tors was made long ago, and the discovery of ferromagnetic 
superconductors (FMSCs) 2,3 displaying coexisting supercon- 
ductivity (SC) and ferromagnetism (FM) well below the Curie 
temperature, has produced a realization of a physically rich 
system that exhibits simultaneously broken 0(3) and £7(1) 
symmetries. Spontaneous breaking of symmetry is respon- 
sible for a wide range of physical phenomena, such as the 
mass differences of elementary particles, phase transitions in 
condensed-matter systems, and even emergent phenomena in 
biology 4 . In many cases, the phenomena can in fact be de- 
scribed by the same equations. Thus, a study of ferromagnetic 
superconductors is of interest not only in terms of having an 
obvious potential for leading to novel devices in for instance 
nanotechnology and spintronics, but also from a fundamental 
physics point of view. 

A successful model describing a FMSC demands that two 
important issues are adressed: i) the symmetry of the pairing 
state, and if) whether the superconducting and ferromagnetic 
order parameters are coexistent or phase-separated in space- 
time. Cooper pairs in conventional superconductors (s-wave) 
are spin-singlets. Thus, s-wave pairing and uniform FM are 
antagonists 5,6 . Spin-triplet Cooper pairs, however, can carry 
a net magnetic moment. Such Cooper pairs are presently the 
prime candidate for explaining the coexistence of FM and SC 
in e.g. UGe2, and URhGe 2,3 . In these materials, SC occurs 
deep within the ferromagnetic phase. It is therefore natu- 
ral to view the SC pairing as originating with electrons that 
also contribute to FM. Thus, the electrons responsible for FM 



below the Curie temperature Tm condense into Cooper pairs 
with magnetic moments aligned along the magnetization be- 
low the critical temperature T c . While spin-singlet supercon- 
ductivity coexisting with uniform ferromagnetism appears to 
be unlikely, it could coexist with helically ordered magnets. 
Tunneling phenomena in such systems have indeed been con- 
sidered theoretically 7,8 . This is, however, physically quite dif- 
ferent from the case we will study in this paper. 

Bulk FMSCs are expected to display an unusual feature, 
namely the spontaneous formation of an Abrikosov vortex 
lattice 9 . Uniform superconducting phases have also been 
investigated 10 , but in a bulk system it seems reasonable to 
assume that this must be associated with a nonuniform mag- 
netic stated 8 ". A key variable determining whether a vortex 
lattice appears or not seems to be the strength of the inter- 
nal magnetization m 11 . The current experimental data 3 avail- 
able for URhGe apparently do not settle this issue unambigu- 
ously, while uniform coexistence of FM and SC appears to 
have been experimentally verified 12 in UGe2. Moreover, a 
bulk Meissner state in the FMSC RuSr 2 GdCu 2 8 has been 
reported^. No consensus has so far been reached concerning 
the correct pairing symmetry for the FMSCs, although theo- 
retical considerations strongly suggest that a non-unitary state 
is favored 1 4 ' 15,16 . In terms of the dt-vector formalism (see 
below), this means that the relation i(dj c x d£) ^ is satis- 
fied, which is equivalent to saying that the Cooper pairs carry 
a net spin 17 . The study of pairing symmetries in unconven- 
tional superconductors has a long tradition 18 where tunneling 
currents have played a crucial role. For the case of spin-triplet 
non-magnetic superconductors, theoretical studies of tunnel- 
ing to a normal metal or ferromagnet have been suggestive 
in terms of establishing the correct pairing symmetry for the 
superconductor?^ .. 

In this paper, we consider quantum transport between two 
thin films of a non-unitary FMSC and an easy-axis ferromag- 
net, respectively. We demonstrate how the resulting conduc- 
tance spectra can be exploited to obtain useful information 
about the superconducting gaps. The SC and FM order pa- 
rameters are assumed not to be phase-separated. Moreover, 
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the choice of a thin film FMSC is motivated by the fact that the 
pair-breaking orbital effect on Cooper pairs with an in-plane 
magnetization is suppressed, and that one will avoid vortices 
present in the compound if the thickness t of the film is smaller 
than the diameter of a vortex 21 , t < £o <C Ao. Here, £o is the 
coherence length (typical size of Cooper pairs) while Ao is the 
penetration depth (typical radius of vortex). In the cases^ 3 - of 
UGe2 and URhGe, this amounts to t ~ 10 nm which is well 
within reach of current experimental techniques. 

We organize our work as follows. In Sec. HJ we establish 
the model and formalism which we will apply to the problem. 
Results are given in Sec. [TTTJ in addition to a discussion of our 
findings. Concluding remarks are given in Sec. [IV] 



II. MODEL AND FORMULATION 

Our model is illustrated in Fig. Q] where the superconduct- 
ing pairing symmetry is equivalent to that of an /12-phase in 
3 He [see Ref . [l7l and the discussion below Eq. (fl}] . 
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FIG. 1: (Color online). Model system of a FM-FMSC interface in a 
slab-geometry in the clean limit. Retroreflection symmetry is not 
broken, since the hole carries the same spin as the incident elec- 
tron. We have sketched gaps corresponding to the analog 1 - of the 
y!2-phase in liquid 3 He. 

The Bogoliubov-de Gennes (BdG) equations for a FMSG 2 ^ 
may be written in compact matrix form 



Mk A k 

A* -M k 



kcr 



Vkcr 



(1) 



where we have introduced the 2x2 matrices 

Mk = fkl - a z U R , 
A k = (A 1 + <x • d k )ia y , 



(2) 



in addition to u k(T = [<i kT<T u k | CT ] T , w k(T = [u kT(T v k j, cr ] T . 
Here, e k is a kinetic energy term, a = (tr x , a y , a z ) is a vector 
containing the Pauli matrices, ?/r,(l) describes the magnetic 
exchange energy of right (left) part of the system. Moreover, 
the d k vectors are given by 



and Ao and A kct/ 3 are the superconducting spin-singlet and 
spin-triplet order parameters, respectively. Finally, it k(T and 
Vka are the wave-function solutions with eigenvalue i? k(T . In 
the following, we set the k-independent singlet amplitude Ao 
to zero, as we do not consider the case of coexistent s-wave 
SC and FM 16 . Also, the opposite-spin triplet pairing giving 
rise to the gap Akfj. is m general believed to be suppressed, 
since it requires interband pairing of Zeeman-split electrons 3 . 
Hence, our model is that of a non-unitary spin-triplet state 
with equal-spin pairing, i.e. Ak-fj. = 0, A kcrcr ^ 0, equivalent 
to the ^42-phase in liquid 3 He (see e.g. Ref. 17) with a non- 
vanishing magnetic moment associated with the Cooper pairs. 

As indicated in Fig. Q] the reflected and transmitted 
electron- and hole-like excitations feel different pairing poten- 
tials due to the orbital symmetry of the superconducting gaps. 
The angle into which they are scattered depends on the spin 
a of the incident electron, since there is a magnetic exchange 
energy Ur present in the superconducting state. In the fol- 
lowing, we will consider the momentum of the quasiparticles 
as fixed on the Fermi surface for spin a so that the supercon- 
ducting gaps correspondingly only depend on the direction of 
momentum, i.e. A kcrcr — > A CT (0J) where 9° is defined in Fig. 
Q] Throughout this paper, we shall insert the superconducting 
gap a priori instead of solving it self-consistently, in order to 
obtain analytical formulas. This is justified by the fact that de- 
tailed calculations taking into account the modification of the 
pair potential near the barrier demonstrate that no new quali- 
tative features appear in the conductance spectra compared to 
when the gap is modelled by a step function at the interface 23 . 
However, the proximity effect at a FM/SC interface may cause 
two important phenomena to occur: i) induction of an SC or- 
der parameter (possible mixture of singlet and triplet) in the 
FM region 24 and ii) the formation of midgap surface states on 
the interface 25 , leading to a suppression of the OP in the vicin- 
ity of the interface. The competition and coexistence of these 
two phenomena has been studed in Ref. 26. In this work, we 
will mainly be concerned with a SC pairing symmetry anal- 
ogous to the A2-phase in liquid 3 He, for which the latter of 
these effects is only present for a specific trajectory of the in- 
coming electrons. Thus, we believe our results capture the 
essential qualitative features of the conductance spectra, al- 
though a self-consistent approach would be required in order 
to obtain the entire picture. 

For the simplest model that illustrates the physics, we have 
chosen a two-dimensional FM-FMSC junction with a barrier 
modelled by the potential V(r) = VqS(z) and superconduct- 
ing gaps A ff (0J,r) = A CT (6»J)e(z). Here, 8{z) and G(z) 
represent the Kronecker delta- and Heaviside functions, re- 
spectively. Solving the BdG-equations and applying the BTK 
formalism 27 , one finds that our system in Fig. Q]is described 
by the wave-functions for particles and holes with spin a in 
the ferromagnet (4> <7 ) and FMSC ( v I' cr ) side of the barrier. 
Explicitly, the total wave-function ^/^ Qt (z) = 0(— z)i/) CT (z) + 
Q(z)^ (T {z) then has the components 
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with 7ff (0) = A a (6)/\A a (6)l 91 + = 9°, Q%_ = 7T - 0J. 
The wave-vectors read fc ' = [2m(i?F + pUl)] 1 / 2 , q a = 
[2m(E F +aUn,)} 1/2 while the spin-generalized coherence fac- 
tors are 

u a {e: ± )K{e° ± )] = \{i + [-]^i-{\K{o° s ±)\/Eyy/ 2 . 

(5) 

In writing down Eq. (@), we have implicitly incorporated con- 
servation of group velocity and conservation of momentum 
parallell to the barrier, i.e. k a sin 9 — q a sin 9°. 

III. RESULTS AND DISCUSSION 

The normal- and Andreev-reflection coefficients can be cal- 
culated by making use of the boundary conditions 

ijf(o) = r(o), 

ii) 2mV r(0) = d/dz[^(z) - ^(z)]\ z=0 . (6) 
Let us introduce Z — 2mVo / fcp and 

r a ± =q° cos6l±k a co$6±ik F Z, (7) 

while Pa (R) = (E F + aU L{R) )/2E F denotes the spin polar- 
ization on the left (right) side of the junction. Our calcula- 
tions then lead to the explicit expressions for the Andreev- 
and normal-reflection coefficients for this FM-FMSC junc- 
tion, namely 

r a e = -1 + 2k a cos 9[u a (91 + )u a {9 a s _){T a + y 

+ v a (e°_ )v a (^ + ) 7 , ip°_ H {e° + ) (Ti )*] ID" , 

rl =4k° cos 9q° cos 8: v49: + )u a (9:_)rM + )/ Da > ( 8 ) 
upon defining the quantity 

D° =u a {9: + )u a {9:_)\T° + \ 2 

- «„(9?-K(^)7.(t)7:(^) x l?-| 2 - (9) 

In the limit of weak FM where all quasiparticle momenta 
equal kp, our results are found to be consistent with a spin- 
generalized augmentation of the equations in Ref. |29l In then- 
general form, the above equations are novel results that in- 
clude the effect of magnetism on the superconducting side 
of the barrier. As demanded by consistency, one obtains to- 
tal reflection \r a e | 2 = 1 when 8 > arcsin(q cr /fc' T ). Having 
obtained the above quantities, one may calculate the conduc- 
tance G(E) of the setup (in units of the normal conductance, 
i.e. A CT (0) — > 0). We find that it is given by 

G(E) =Y J G a (E)/Y J F a (10) 



where we have defined the quantities 

,tt/2 

G a (E)[F a } = / d9 cos 9G a (E,9)[f a }P^P^, 

J-n/2 

G a (E,0) = 1 + K(E,9)\ 2 - \r°(E,9)\ 2 , (11) 

and f a = 1 - |1 - 2fc CT cos6>/Y^| 2 . We next demonstrate 
how the conductance spectra yields useful information about 
the superconducting gaps upon varying the relative orientation 
of the gaps, their magnitude, and the strength of the magnetic 
exchange energy on each side of the junction. To be specific, 
we first first consider the analog of the A2-phase in liquid 3 He, 
such that the gaps may be written 2 ^ 

A T (0J ± ) = -A T ,oe^±-«), A^) = A lfi e-^-?\ 

(12) 

as illustrated in Fig. Q] We stress that 6° ± is not the global 
broken U(l) superconducting phase, but rather an internal 
phase originating with the odd symmetry of the p-wave gaps 
(see Fig. Q]). The exchange energy of the FMSC will be kept 
fixed at [/r = 0.05Ef- In Fig. [2] we have plotted the con- 
ductance spectra for a FM/FMSC junction for three distinct 
cases. We have defined the ratio between magnitude of the 
superconducting gaps as Ra — Af.o/Aj.o- Introducing the 
dimensionless barrier strength Z = 2mVo/fcF, where fcp is 
the Fermi momentum, we consider i) the metallic contact case 
with no barrier (Z = 0), ii) the intermediate case with a mod- 
erate barrier (Z = 3), and Hi) the tunneling limit with an in- 
sulator in the junction (Z — » oo). For each case i)-iii), we 
have allowed the magnetic exchange energy £/l to vary from 
weak, favoring f-spins (U — Q.05Ep) to strong, favoring f- 
spins in one case (U — 0.5Ep) and j-spins in the other case 
(U = — 0.5-Ef)- These are shown in descending order in each 
column. We have also considered the conductance G(E) for 
several values of {a, (3}. For the gaps chosen above, we find 
that the G{E) did not depend on different choices of these 
parameters. This can be understood by observing that the an- 
gular averaging over G IJ {E,9) allows for factors e" 1 "^ to 
be separated out on equal footing as the factor corresponding 
to the globally broken U(l) symmetry, whose value does not 
affect the conductance spectra. This will, however, not be the 
case for other possible triplet symmetries, as discussed below. 

An important, and obvious, feature of the conductance 
spectra is a characteristic behavior occuring at voltages cor- 
responding to E = A a fl, a =|, I, where peaks are displayed. 
This offers the opportunity to utilise the conductance spectra 
to reap explicit information on the size of the superconduct- 
ing gaps in the FMSC. From Fig. |2j it is seen that the ef- 
fect of increasing the exchange energy on the ferromagnetic 
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side to J7l = ±0.5£f is a sharpening of the peaks located at 
E = A a o, where a is the spin-species energetically favored 
by [/l- Concomitantly, the peak located at E = A_ CT .o is sup- 
pressed. Such a response is consistent with what one would 
expect, since increased/decreased spin polarization on the fer- 
romagnetic side enhances/suppresses the conductance of the 
corresponding spin component. In the tunneling limit (large 



Z), we see that the conductance spectra exhibits sharp transi- 
tions at E = A CT) o, corresponding to the sudden appearance 
of a tunneling current as the voltage exceeds the magnitude of 
the gaps. We have also considered G(E) in the case of vanish- 
ing FM on the left side, i.e. a N/FMSC junction. The results 
we find are very similar to the case of weak FM displayed in 
the upper row of Fig. |2j and are therefore not displayed. 



i) Metallic contact (Z = 0) 



'.) Intermediate case (Z — 3) 





Hi) Tunneling limit (Z — > co) 
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FIG. 2: (Color online). Plot of the conductance G(E) for a FM/FMSC junction in the case of a metallic contact, the tunneling limit, and an 
intermediate case. The gaps are given by Eq. J 1 2b in this case, for which G(E) is found to be insensitive to {a, f3}. In the above, {a, /3} = 0. 
The columns for each case provide the spectrum for Ul = 0.05-Ef, Ul = 0.5-Ef, and Ul ~ — 0.5Ef, in descending order. For each figure, 
we have plotted several ratios between the magnitude of the superconducting gaps. These are given by Ra = {2, 3, 4, 5}, represented by the 
magenta, blue, green, and red line, respectively. 



For the gaps in Eq. ( TTSI i. the conductance was found to 
be insensitive to {a,/3}. However, in general this is differ- 
ent, and the orientation of the gaps relative to the barrier is 
essential in determining the conductance spectrum. We il- 
lustrate this with a somewhat different choice of anisotropic 
gaps. When the superconducting gap is oriented in a fash- 
ion that leaves it invariant under inversion of the component 
of momentum perpendicular to the junction, k z — > — k z in 
this case (equivalently 9 — > it — 9), no zero-bias conductance 
peak (ZBCP) should be expected as there is no formation of 
current-carrying zero energy states 25 . However, when the gap 
changes sign under such an inversion of momentum, ZBCPs 
appear 29 . Since the momentum perpendicular to the junction 
of the hole-like excitation in the anisotropic superconductor is 
reversed compared to the electron-like excitation, a gap that 
satisfies A a (9) ^ A a (n — 9) will cause the hole to feel a dif- 
ferent pairing potential than the electron-like excitation. This 
is the fundamental reason for the appearance of a ZBCP. How- 
ever, in the present case of p-wave superconducting gaps anal- 
ogous to the A2-phase in 3 He, both gaps are asymmetric un- 



der the operation 9 — > it — 9 regardless of the value of {a, f3} 
except for the single value 9 = 0. Therefore, a small contri- 
bution to zero energy current-carrying states, i.e. G(0) ^ 0, 
will occur as shown in Fig. [2] However, this contribution does 
not lead to a ZBCP, for which all 9 (see Fig. |2]i contribute to 
the formation of zero energy states due to the asymmetry of 
the gaps. In a model where the p-wave gaps are represented 
by the odd-symmetry analog of d x 2_ y 2-wave gaps, i.e. 

A T (0) = A T , cos(0 - a), A;(0) = A;, cos(9 - (3), (13) 

the formation of ZBCP will then be predictable from the 
orientation of the gaps as these can now display symme- 
try/antisymmetry/asymmetry when 9 — ► tt — 9 7 depending on 
{a, (3} (see Fig. [TJ. Indeed, insertion of the above gaps into 
our model produces conductance spectra that display a ZBCP 
for e.g. a — (3 = 0, as can be seen in Fig. [3] In this case, as 
in Fig. |2] the conductance spectra also allows for the magni- 
tude of the superconducting gaps to be read out, although the 
features are not as clear as those seen in Fig. [2] 
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FIG. 3: (Color online). Plot of the conductance G(E) for FM/FMSC 
junction with gaps given by Eq. d 1 3h . for which G(E) is sensitive 
to {a, f3}. We have chosen Z = 3, and Ra = {2, 3, 4, 5} are 
represented by the magenta, blue, green, and red line, respectively. 

From the results of Figs. [2] and [3] it is clear that the con- 
ductance spectra G(E) exhibit strong dependence on the ex- 
change energy, while the relative orientation of the gaps is 
insignificant for the superconducting gaps Eq. (fT2] i. Thus, our 
results will shed light on the magnitude of the various com- 
ponents of the superconducting gaps and their relative orien- 
tations in k-space if the gaps display symmetry/antisymmetry 
under 6 — > it — 9 for some orientation. Moreover, we are 
dealing with an easily observable effect, since distinguishing 
between the peaks occuring for various values of Ra requires 
a resolution of order 0(lO -1 A^o), which typically corre- 
sponds to 0.1 — 1 meV. These structures should be then readily 
resolved with present-day STM technology. However, in or- 
der to do justice to the experimentalist, it should be pointed 
out that a challenge with respect to tunneling junctions is 
dealing with non-idealities at the interface which may affect 
the conductance spectrum. Also, the importance of spin-flip 
processes in the vicinity of the interface (if such are indeed 
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present) has recently been pointed out 30 . 



IV. CONCLUSION 

In summary, we have studied the conductance spectra 
G(E) for systems consisting of a ferromagnetic metal sepa- 
rated from a non-unitary p-wave FMSC by a thin, insulating 
barrier. We have considered the cases of weak, intermediate, 
and strong barriers, and considered three different strengths 
of the ferromagnetic exchange energy. We have considered 
two classes of anisotropic spin-triplet superconducting gaps, 
with results given in Figs. (O and (01. Our results show how 
the magnitude of the superconducting gaps A a , a =f , J, may 
be inferred from the conductance spectra. Moreover, the class 
of superconducting gaps given in Eq. (fT2] > renders G(E) 
insensitive to {a,f3}, the results are shown in Fig. d2j fol- 
ia, (3} — 0. Conversely, the orientations of the p-wave gaps 
modelled by Eq. (1131 , specific values of {a, (3} are essential 
to the formation of ZBCPs in G(E) in Fig. ©, and the 
characteristic behavior at E = A ff o, o = T>I- These results 
should provide useful insights in determining both the relative 
orientation between the superconducting gaps, as well as their 
magnitude, in ferromagnetic spin-triplet superconductors. 
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